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Abstract 

We consider the Ising spin glass for the arbitrary spin S with the short- ranged 

interaction using the Bethe- Peierls approximation previously formulated by 

Serva and Paladin for the same system but limited to S = 1/2. Results 

obtained by us for arbitrary S are not a simple generalization of those for 

5 = 1/2. In this paper we mainly concentrate our studies on the calcutation 

of the citical temperature and the linear susceptibility in the paramagnetic 

phase as functions of the dimension of the system and spin number S. These 

dependences are illustrated by corresponding plots. 
75.10.Nr, 75.50.Lk 
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I. INTRODUCTION 



The study of glasses is today one of the most relevant and actual problem in condensed 
matter physics. In particular research around spin glasses in finite dimensions is very active 
since it is still unclear if they share some the qualitative features of the mean- field theory 
of the model Sherrington- Kirkpatrick (SK) 0^ However, there are recent studies 
which indicate difficulties to extend the molecular field approximation (MFA) scenario to 
realistic spin glasses with short- range interaction and decide "a priori" which properties 
survive and which must be appropriately modified. Recently, in an interesting paper M , an 
approach beyond the MFA has been achieved for an d- dimensional Ising SG model (S = 1/2) 
with short- range interactions on a real lattice using an extension of the Bethe- Peierls 
approximation (BPA) |7| to the spin glass problem via the replica trick. This approach 
seems to be very promising to estabilish a direct contact with the results obtained by different 
authors for the infinite- ranged version and to controll possible deviations for short- ranged 
glasses from the well acquired MFA scenario. Quite recently the Parisi ansatz has 
been investigated for the Ising spin glass with S = 1/2 using the generalized of the Bethe- 
Peierls metod named by the authors "a variational approach" where finite clusters of spins 
interacts and the sample averaging is properly taken into account. The result is qualitatively 
similar to that obtained in the frame of the MFA with some quantitatively modifications 
due to short- range order interactions. 

All studies mentioned above has been performed for the standard Ising model with 
S = 1/2. Therefore it seems to be quite interesting to extend the methods applied for the 
Ising spin glass, where the number of spin is arbitrary. In the spite of numerous works on 
the Ising model S — 1/2 a little attention has been devoted to the same system with S 
arbitrary. In the spin glasses solutions for the arbitrary S is not a simple generalization of 
that for the Ising model with S = 1/2. The reason is that for higher spins (5$)™ ^ Si or 
const which leads to parameters diagonal in replica indices. This considerable affects the 
results for the S = 1/2 Ising SG. 
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We will use the Bethe- Peierls approximation (BPA) with some necessary modifica- 
tions. As we will see the critical temperature depends as well as of the dimension of the 
system and of the spin number. In addition, linear susceptibility even in the paramagnetic 
phase are nontrivial functions of the temperature. 

The paper is planned as follows: In Section [TT] we shortly present the idea of the Bethe- 
Peierls approach for the the Ising spin glass @ and give equations for the critical temperature 
and the linear susceptibility in the paramagnetic phase when the spin number S is arbitrary. 
The Section |TTl| contains plots as the critical temperature depends of the dimension and the 
value of the spin number. Finally in Section |IV| some conclusion are drawn. 



II. BETHE- PEIERLS METHOD FOR SPIN GLASSES 

The our starting point is the Hamiltonian of the Ising model with the arbitrary spin: 

H = —-Y^ijJijSiSj , (1) 

where Si = Sf is the z-component of the spin. As usual eigenvalues of Si runs from —S to 
S, where S is arbitrary. In (|l|) the summation over i,j comprises only nearest neighbours. 
Interaction parametrs Jjj are variables obeying random distribution that is , dichotomic, 
gausssian, etc. For the simplicity, in order do not complicated the main idea, we will assume 
that Ji j is a random variable with Jjj = ±J with the equal probability for + and — sign. 
Using the replica trick the free energy can be written as follows: 

(2) 

where [• • -} av denotes a sample averaging and H a is the replicated Hamiltonian (JTJ) . 

Working directly on the real lattice, the basic idea of the BPA for spin glasses |7j is 
to take into account the correct interactions inside replicated clusters (cl), constitued by 
a central spin S and its 2d nearest- neighbours Sf, {i = l,---,2d}, and to describe the 
interactions of the cluster borders with the remnant of the system by means of effective 



f3F = lim^-Trexp -^VH £ 
n — ' 
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couplings among replicas to be determined self - consistently. Therefore the free energy in 
this approximation can be written as follows: 



- f3F = lim n ^ - In Z n 

n 



(3) 



where 



Z n = K((3)Tr 

i?J 2 2,1 



2<l 



i=l.a 



/ ^2 j2 2d n 

6XP \~2~ 



(4) 



x exp 



2d \ | 



In Eq. 



a = 1 • • • n denotes the replica indices, K (f3) is a multiplicative constant which 
depends on the temperature but not on the lateral spins of the cluster, [i a ^ a i and /i = // a , a , 
according to Bethe- Peierls ansatz, describe the interaction betwen the "external world" 
and the lateral spins of the replicated cluster. The difference between Eq. (Q) and the 
corresponding formula of Ref. J7[ is that in (|J) we have the additional parameter /x. This is 
a consequence of the fact that for an arbitrary spin S\ a ^ 1/4. 

Effective couplings fi aa i and \i are calculated from the following equations: 



{Si ta Si ta ') — (So,aS(),a') 



(5) 



with i — 1 • • • 2d, where 



(...) 



Trexp {-(3H eS ) ■ ■ 
Trexp (-(3H eS ) 



(6) 



with 

H e f[ 



In 



exp 1/3 Jo,iSo, a Si, a 



2d 



i,a aj^a' i=l 

Above and at the critical point the spin glass parameters (Sk, a Sk,a') = Q a ,a' — for 
a/a' (but not for a = a'). It is easy to calculate that the effective couplings fi a ^a' obey 
the same conditions. If we are interested in the calculation of the critical temperature it is 
sufficient to formulate the equation (^j) to the lowest order in ^i a -^ aP rime . After straighforward 
algebra we get for a ^ a' the following result: 



2d 



(Si,aSi,a') ~ P 2 J 2 f^a,a' X] [(^^j')o] 

3=1 



(8) 



with i, j — 1 • • • 2d and 



2f/ 



I av ' 



(9) 



i=l 



where a ^ a' and 



(•••)c 



Tr exp (p Egi Jo.i-S'o-Si + E -=i S? 



Tr exp (/? Egi Jo^o^ + ^/i E^i Sf 



(10) 



/-^n^i^Trexp 


'/?E-=i (j ^ + V / V)^ 




roo^l^Trexp 







with 



T>Xi = — exp (—x 2 /2) dxi 

V 27T V 7 



(11) 



Thus due to the translational symmetry for averaged over disorder correlation functions the 
equation for the critical temperature takes the form: 



(St)t „ +(2d-l) (SiSM 



2d 



(SoSi)l 



(12) 



where i 7^ j numbers of arbitrary lateral spins of the cluster. Additionally we must to take 
into account the equation: 



(13) 



With Eqs (|9]) and (pT0| ) after detailed calculations we can formulate the equation for the criti- 
cal temperature in terms of functions F to F 4 which depend on the temperature, dimension 
of the system, and the spin number. The sample averaged correlation function with i 7^ j 
has the form 



(SiSj)l 



n 



av Fi ' 



(14) 



whereas 
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av F j 



3_ 

2 ' 



F 2 



av F 2 



and 



F 2 



av F $ 



The form of funcions Fi (i = • • ■ 4) is following 
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M=-5 



PxQ s [/? (M + /i 1/2 



X 



2d 



F 1 = S 2 E / VxB 



x 



M=-5 
oo 



(M + /i 1/2 x)] Q 5 [/3JS (M + /i 



,V 2 



£>xQ s [/? J5 (m + ^ 

oo L ^ 



V2, 



2d-2 



F 2 = S E ( /°° ^{^s [/3</S(M + /i 1/2 x)] +5| [/3JS(M + // 1/2 a;)]} 



M=-S 



Qs [PJS (M + fi 



1/2 



Pa;Qs [PJS (M + /i 



1/2, 



2d-l 



I ( poo 

F 3 = S E M j/ [/3J5(M + /i 1/2 a;)] g s [/3J5(M + ^ 



M=-S 



,1/2 



;r 



oo 



£>xQ s f3JS[M + n 1/2 



x 



2d— 1 



and 



^4= E M2 

Af=-S 



PxQs [/3 JS (M + // 



1/2 



;r 



2J 



In Eqs 



Qs (y) 



sinh 




(1+ 
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sinh (§) 





6 



B s (y) 



2S Ul,th 



is the Brillouin function, and 



B' s {y) 



1 + 2S^ y 



dB' s (y) 



(24) 



dy 



(25) 



Taking into account Eqs fll2|) and (|T3| ) together with (|14|)- (0) we can write equations for 
the critical temperature as follows: 



(2d - 1) F? + F* = 2dF* 



(26) 



and 



Fo — Fa 



(27) 



Obviously the solution of (|26| ) and ( p7[ ) needs numerical calculations. 

Which concerns the linear susceptibility in zero magnetic field we define it as follows: 



X 



N 



N 

E 



d(Si)T,h 
dh 



\h=0 



(28) 



where (• • ')t,h denotes the thermal averaging with the Hamiltonian ([I]), when the term 
—h St is added. A first step is to calculate the susceptibility in the paramagnetic phase 
with local magnetizations (Si)T = 0, where (Si)T = (Si)r,h=o- In that case 



N 

x = |fE [ W>t] 

i,3 



It is easy to show that for h = and the symmetric probability distribution for J, 



'•j 



[(SiSj) 



av 1 >3 



After some calculation (see, for example Ref. |Kj) we get that 



X = P (Sl)o 



(29) 



(30) 



(31) 



with k = • • • 2d. 



III. RESULTS 



Our results are illustrated by plots [T]- |j. In Fig. [1] the dependence of the citical tem- 
perature T c (in units of the constant J) scaled by \/2d of the dimension of the system d for 
a few values of spin S = 1/2, 1, 3/2, 2 and 3 are given. The larger the spin, the higher the 
corresponding line. In Fig. Q variations of the T c j\[2d with spin number for d — 2 ( the 
lower line) and d = 3 (the upper line) are plotted. It is seen that the critical tempetature, 
in general, increases with the increasing of the spin number but this dependendence cannot 
be represented in an explicit form and is more nontrivial comparing to the simple magnetic 
systems as, for example, for a ferromagnet where T c ~ S. The reason about scaling of the 
ctitical temperature by \pld~ will be explained in the next Section. On Figs ||| and § we show 
the dependence of the linear susceptibility in the paramagnetic phase of the temperature 
(J = 1) for d = 2 and d = 3, respectively. The values of S are 1, 3/2 and 2. The larger the 
spin the higher the corresponding line. Obviously the lines in Figs |3] and |] terminate at the 
critical temperature since to calculate x below T c we must enter into theory the spin glass 
order parameters. At present our purpose is only to show that the linear susceptibility, even 
in the paramagnetic phase, is a nontrivial function of the temperature for spins higher than 
S = 1/2. From plots |3] and || it is seen a tendency that x increases with increasing of S. 



IV. FINAL REMARKS 

The studies of the spin glasses with the short ranged interaction is undoubtetly a difficult 
problem among the theories of amorphous systems since the complicated nature of the 
randomness interplays spatial correlations of spins. At the present practically there is no 
developed a systematic method to investigate such systems, as in the case of the long ranged 
(more strictly infinite- ranged) Sherrington- Kirpatrick type models where MFA is valid 
|H||]. It is expected that the BPA is able to give more aacurate estimation of the critical 
temperature for the spin glass systems with short range iteraction than the MFA. Which 
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concerns our problem a natural question arises which is result when the dimension of the 
system is ifinite. It can be easy shown that if d — *■ oo one obtains the Sherrington- Kirpatrick 
theory for the Ising spin glass with an arbitrary spin. After rescaling Jjj — > Jij/\r2d, 
J — > J/V2d and changing /z Qj< y = 2dq a ^ a ', p> = 2dp with p = q aj0l proceeding in the similar 
line as in Ref. J7| one obtains 

q a ^ a ' = (S a S a >) (32) 

and 

P = (Si) , (33) 

where 

_ Tr exp (gg q«, a >S a S al + pE a g) ; • • , /2 x 

Tr exp ^^g- Eq^q' q a ,a'S a S a ' + P Ea ^« J 

where S a is the 2- component of the spin operator referred to an arbitrary site. Hence 
in Figs [l] and |2| in order to obtain resonable results for higher spins we scaled the critical 
temperature, expressed in units J, by V2d. We are aware that our consideration is a first 
step to recognize some properties of the Ising SG with an arbitrary spin. It would be 
interesting to obtain the properties of the system in the SG phase at least in the replica 
symmetric theory. This is a complicated task even for S — 1/2. Therefore, further work will 
be necessary to elucidate this problem. 
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FIGURES 

FIG. 1. The variation of the critical temperature T rescaled by the factor (2d)~ 1 ^ 2 with the 
dimension d of the system for spin numbers S = 1/2,1,3/2,2 and 3 . The larger the spin the 
higher the corresponding line. Here J = 1. 

FIG. 2. Rescaling critical temperature T c / (2d)~ 1 ^ 2 versus the spin S for d = 2 and d = 3 
marked by lower and upper lines, respectively. Here J = 1. 

FIG. 3. Linear susceptibility % i n the paramagnetic state as a function of the temperature in 
units of J for d = 2 and 5 = 1, 3/2 and 2. The larger the spin the higher corresponding line. 

FIG. 4. Linear susceptibility \ i n the paramagnetic state as a function of the temperature in 
units of J for d = 3 and 5=1, 3/2 and 2. The larger the spin the higher corresponding line. 
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